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Abstract 

We show that many of the unusual properties of the one-dimensional random 
quantum Ising model are shared also by dilute quantum Ising systems in the 
vicinity of a certain quantum transition in any dimension d > 1. Thus, 
while these properties are not an artifact of d = 1, they do require special 
circumstances in higher dimensions. 
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There is considerable current interest in the properties of phase transitions in random 
quantum systems. Experimentally accessible quantum transitions, such as the transition 
from an insulator to a metal or a superconductor [Q, and various magnetic-nonmagnetic 
transitions in heavy fermion compounds ||^, high- Tc cuprates 0], or insulating dipolar Ising 
magnets often occur in situations with strong randomness, and are only poorly un- 
derstood. Theoretically many authors have analyzed the random Ising chain in a 
transverse field - perhaps the simplest random quantum system. In particular. Fisher |^ 
used a real space renormalization group (RG) approach to obtain a detailed description of 
the thermodynamics and static correlation functions in the vicinity of the critical point. The 
properties of the system were found to be very unusual as compared to conventional quantum 
critical points. Specifically, length scales were found to diverge logarithmically with energy 
scales, implying a dynamic critical exponent z = oo. The transition was shown to be flanked 
on either side by "Griffiths" regions with a susceptibility diverging due to contributions from 
statistically rare fluctuations. There are very few reliable results on other random quantum 
transitions, especially in finite dimensions d > 1; thus it is important to understand if the 
anomalous properties of the random quantum Ising chain are a specialty of (i = 1, or if there 
exist higher dimensional quantum transitions which share these properties. Numerical work 
on higher dimensional transverse field Ising spin glasses has found evidence for the presence 
of Griffiths regions, but the dynamic scaling at the critical point seems conventional 0. 

In this paper, we provide a simple example of such anomalous scaling in higher dimen- 
sional random quantum Ising systems. We consider bond or site diluted Ising models with 



short range interactions. As was first suggested by Harris [|10|, and as we argue below, these 
models have two quantum transitions (See Fig |l]): At low dilution, below the percolation 
threshold, there is a phase transition when the long range ferromagnetic order is destroyed 
by increasing the transverse field. This is expected to be in the universality class of the 
generic random bond quantum Ising transition. Right at the percolation threshold, there 
is a finite range of transverse field strengths at which the system remains critical. There 
is thus another quantum transition, across the percolation threshold, at low but non-zero 



transverse field strengths, which is potentially in a different universality class. These two 
critical lines meet at a multicritical point (M). The properties of the second transition, at 
the percolation threshold, are determined largely by the statistics and geometry of the perco- 
lating clusters about which much information is available. This permits us to make definitive 
statements about the scaling properties of this transition in any dimension. We show that 
the dynamic scaling is activated, with length scales diverging logarithmically with energy 
scales. We also demonstrate the existence of Griffiths phases on either side with diverging 
susceptibility. Our approach shows clearly the connection between these behaviors, and 
T = phase transitions at which quantum fluctuations are "dangerously irrelevant" []8|,pT|; 
indeed, various exponents of the transition are given by those of the classical percolation 
transition, and are hence known either exactly or numerically. We also obtain bounds on 
exponents characterizing the multicritical point. 

For concreteness, we consider bond-diluted Ising models defined by the Hamiltonian 

^ = -E^^.<^|-E^<' (1) 

where (ij) labels the nearest-neighbor sites of a c?-dimensional lattice {d > 1), and Jij equals 
with probability p and equals J > 0, with probability 1 — p. The transverse field h is 
taken to be non-random, although our results remain valid for weakly random distributions 
as well. At p = 0, a.s h is increased, there is a zero temperature phase transition from 
an ordered ground state to a disordered one (See Fig |l]). On the other axis, when h = 0, 
so that the system is classical, there is a percolation transition at p = Pc- For p < Pc, 
there is a thermodynamically large connected cluster, while for p > pc, there only are finite 
connected clusters. For p < pc, for small enough h, the system retains long range order. 
This is ultimately destroyed for some h > hc{p), with hc{p) expected to be a monotonically 
decreasing function of p. On the other hand, if p > Pc, there is no long range order for any 
h. 

Now consider p = pc- Again there is zero magnetization and no long range order for any 
h. However the system stays critical for h < Hm = hc{pc) (Fig 0). To see this, note that, 
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although there is no thermodynamically large connected cluster at pc, there still is an infinite 
connected cluster with a fractal dimension smaller than the spatial dimension. The spins on 
this cluster align together at h = 0. A small but non-zero h is not sufficient to destroy this 
order on the critical cluster. In fact, two spins on any sufficiently large finite cluster remain 
strongly correlated with each other for small h. The critical cluster eventually loses order 
when h reaches Hm- To make this more precise consider the disorder averaged, equal time (r), 
two point spin correlation function. Spins at points and x are correlated only if they belong 
to the same cluster; thus G{x) = [{a'{x,T = 0)or^(0,0)) - {a'{x,0)) {a'{0,0))] = J dCCV{0 
and X belong to the same cluster) ^(if and x belong to the same cluster, they have 
correlation C), where angular brackets are averages over quantum/thermal fluctuations, 
square brackets represent disorder averages, and V{E) is the probability of the event E. 
At p = Pc, the flrst probability ~ for large x |]12[. The arguments above imply 



that the integral over C is non-zero and independent of x for large x for h < Hm- Thus 
G{x,p = Pc) ~ aj-'i+s-'yp fQj, large x for h < Hm, and there is a critical line for h < Hm at 

P = Ve- 
in, this paper, we will primarily focus attention on the transition across this critical line. 
We show that at T = 0, the properties of this transition are strikingly similar to the d = 1 
results for the random quantum Ising chain First consider the equal-time two point 
function away from pc- Arguing as in the previous paragraph, we get G{x) ~ P(0 and 
X belong to the same cluster) ~ x~'^^'^^'^p f {x / S,) for large x where the correlation length 
^ ~ \p — Pc\^'^^- Similarly, for p < Pc the mean magnetization [(cr^(x,0))] ~ P(any given 
site belongs to inflnite cluster) ~ {pc — pY"^. The exponents Pp,rip, z/p are those of classical 
percolation theory, and are known exactly in d = 2 and for d > 6. For any d < 6, they 
satisfy Pp = Vpi^d — 2 + ?7p)/2. 

Now consider dynamic correlations. The imaginary part of the local dynamic suscepti- 
bility [Xi,(^)] = SArP(given site i belongs to cluster of sites)/ ydyViii site i belongs to 
cluster of A^ sites, 'x'liyj) = y). The energy levels of a cluster of A^ sites can be described for 
/i ^ J as follows: The two lowest levels correspond to the states of a single effective Ising 



spin with magnetic moment ~ in an effective transverse field hef[,N- For large N, hes,N can 
be estimated in A^*'' order perturbation theory to be hexp{—cN). (In general, there would 
also be a prefactor that varies as a power of A^; we will ignore this as it is subdominant 
to the exponential in N. When necessary, we will indicate the modifications induced by 
keeping this prefactor). The quantities h and c are of order h and ln(J//i) respectively but 
their precise values depend on the particular cluster being considered. As the distribution of 
h and c is not expected to become very broad near the transition [|l^, we will replace them 



by their typical values ho and cq respectively. Apart from these two lowest levels, there are 
other levels separated from these by energies ~ J. These can be ignored for the low energy 
physics, and for small u <^ h, we only need to consider large clusters. 

We now need the following results [|1^] from percolation theory for the probability P{N, p) 
that a given site belongs to a large cluster of A^ sites. At p = Pc, P{N,pc) ~ A^^^^ where 
T equals 5/2 for d > 6, and equals 1 + d/D for d < 6 where D is the fractal dimension of 
the critical clusters. In particular r = 187/91 in = 2, and is approximately 2.18 in d = 3. 
Away from pc, for d < 6 or d > 8, P{N,p) satisfies the scaling form 



P{N,p)^N'-^g[N/^^) (2) 

The scaling function g{y) is universal but is different for p < Pc and p > Pc] it approaches 1 
for y -C 1, while for y ^ 1, 

g{y,P>Pe)-^y-'^^e-'^'' 

g{y,p<p,)^y-''+^e-^-y'-'''' (3) 

with 6=1, 3/2 and 5/2 for c/ = 2, 3 and > 8 respectively, and 6' = 5/4, -1/9 in = 2, 3 
respectively. The constants c+,_ are of order unity. For 6 < < 8, P{N,p) satisfies a more 



complicated two- variable scaling form |Tj]. For simplicity, we shall not discuss this case 
here, though including it is straightforward. 

Using these percolation results, we get for the dynamic susceptibility of the Ising model 
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Note that the scahng variable is / u) / S,^ . This is a precise statement of the activated 
dynamic scahng mentioned earher, which has thus been shown to occur in all ci > 1 in 
the present model. Asymptotic forms in various limits can be obtained from the limiting 
behavior of g{x) described above. For p > pc, we get x'l ~ uj^^~^°'{ln{hQ/uj)Y^'^ with a ~ C,^^ 
so that a.t p = pc,a = (Including the power-law prefactor in /iefr,Af will only change the 
power of ln(l/co') in by ~ C,~^, and similarly for the prefactor in the expression for p < Pc 
given below). Note that just on the disordered side, the system is gapless with a power-law 
density of states. The physical origin of this is, as usual, the presence of rare, large clusters 
with arbitrarily small energy gaps. For p < p^, the presence of the infinite cluster (and 
the associated long range order) gives rise to a delta function at = 0. For u ^ 0, x'Li^) 
is still determined by contributions from the finite clusters. Proceeding as before, we find 
X'[{uj ^ 0) ~ (l/cj)(ln(/io/cu))i-"exp (-/t(ln(/io/cj))^"^/'^) with k ~ ^-^(i^iM. Again the 
system is gapless. The gaplessness of both the ordered and disordered phases in the vicinity 
of the transition is unlike quantum transitions in pure systems, but is probably generic to 
many random quantum transitions [P,p!5|,p| . 

The magnetization in response to a uniform external applied magnetic field H along the 
z direction can be calculated similarly. For small H <^ h, only large clusters contribute. 
The magnetization per site of a cluster of size is that of an Ising spin of magnetic moment 
iV in a transverse field hes,N, and is therefore given by 

NH 

""^"^ = ((iVi?)^ + /4,.)"^ 

Thus the total magnetization per site (after subtracting the regular contribution of the 
infinite cluster for p < p^) is 

MiH) - M{H = 0) ~ / dNj^^g{N/e)MM{H) 

The singular part therefore has the scaling form 



with c a non- universal constant, and $(?/) a universal function which is related to g{y) by 
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Again note that the scaling variable is \n{hQ/H)^~^, which is unlike conventional critical 
behavior, but is similar to the d = 1 random quantum Ising chain 
We thus have the following asymptotic forms as if 0: 

(In(V^))'"" P = Pc 

e{Hho/H)y-%H/hoy p>pc (6) 

^D{i-i/d)Qj^(/^^/^))-e'+i-i/dg-{yin(ho/H))i-Vd p < 

with 7,7' ~ (as with the dynamic susceptibility, including the power-law prefactor in 
hcs^N changes the power of ln(l/ii) for p > pc, and the prefactor for p < pc by ~ C,'^)- Note 
that the magnetization rises as a power of H, with a continuously varying exponent which 
is smaller than 1 in a region of the disordered phase close to the transition. This is again 
similar to the d = 1 result, and gives rise to a divergent linear susceptibility throughout 
this region. In the ordered side dM/dH ~ ^/H with weak corrections. Thus the linear 
susceptibility diverges in the ordered side as well. 

A further similarity with the d = 1 results can be found by generalizing the problem to 
systems with Potts symmetry. For the random g-state quantum Potts model in c? = 1, it has 



been shown [|T6| that all the critical exponents and the scaling functions for the mean spatial 
correlations and magnetization are independent of q, i.e, they are the same as for the Ising 
case. For the diluted models that we have considered here, it should be clear that a vertical 
critical line exists at the percolation threshold for all q. The exponents and appropriate 
scaling functions of the corresponding quantum transitions are again independent of q as 
they are determined mainly by the geometric properties of the lattice near the percolation 
threshold. As in d = 1, all the q dependence is in non-universal quantities and in a high- 
energy cutoff limiting the regime of universal scaling behavior. 



So far, we have focussed on the T = scahng properties, and demonstrated the similar- 
ities with the d = 1 results. What about the finite T properties ? For the classical dilute 
Ising model ai p = pc, the correlation length at finite T behaves as exp(contant/T) ||10|| . 
This is essentially due to the presence of one dimensional segments in the critical percolating 
clusters. For the quantum problem ioi h < Iim, these one dimensional segments would give 
rise to a thermal correlation length (^t) with a similar exponential dependence on 1/T, and 
a prefactor that is a power-law in T; this is the behavior in the "high T" region of Fig |^. 
In contrast, for the d = 1 random quantum Ising chain, the correlation length rises as a 
power of ln(l/T) at the critical couplings. Away from the critical point, the crossovers are 
as shown in Fig|^. The low T behavior appears when ^ ~ ^t, or T ~ \n~^ (1 / (\p — pd) ■ On 
the disordered side, the low T system is described well as a collection of rigid Ising clusters 
with effective transverse fields and a size distribution as before; this leads, for instance, to 
a linear susceptibility xt ~ T~^~^'^ (up to log corrections) with k, ~ ■ On the ordered 
side, there is a finite temperature phase transition; as in the classical case, as p Pc, the 
transition temperature falls to zero as ~ ln~^(l/(pc ~ p))- 

We now turn to the special point M (Fig |I|). Correlations at the point M should decay 
faster than along the vertical critical line considered above. Thus, if Gm{x) ~ x~'^, then 
(f) > d — 2 + T]p. Similarly, as p approaches pc from below at fiM, the magnetization should 
go to zero faster than for h < Iim- Thus Pm > Pp- A detailed theory of the scaling at this 
point we leave as an open question. 

Finally it is interesting to ask if this vertical critical line at the percolation threshold 
continues to exist for systems with continuous symmetry, such as the 0{N) rotor models. 
For N > 2, the presence of one dimensional segments in the critical clusters implies rapidly 
decaying correlations in large clusters. Thus the vertical critical line will not be present for 
0{N) models for N > 2. The 0(2) case is special as correlations only decay as a power-law 
in one dimension. Whether this is like the Ising or the 0{N > 2) system we leave as another 
open question. 

To summarize, we have presented a simple example of a random quantum transition 

8 



in dimensions d > 1 which exhibits many of the properties of the transition in the d = 
1 random quantum Ising chain. In particular, the dynamic scahng was activated, with 
ln(l/energy scale) ~ C,^. Further, there were Griffiths regions on either side of the transition, 
with a singular density of states and a diverging susceptibility. It would be interesting to 
find experimental systems where this transition can be studied. (One possibility is the 
system LiHoxYi-x^A where, however, the presence of dipolar interactions may complicate 
the situation.) Theoretically, an important feature of this transition, as in the d = 1 system, 
is that it is controlled by a classical fixed point with quantum fluctuations being "dangerously 
irrelevant". This feature is also found in quantum Ising models in a random longitudinal 
field, which undergo, for d > 2, a. phase transition from an ordered to a paramagnetic phase 



at T = [T^. In fact, this fixed point also controls the non-zero T phase transition, where it 
has been argued that the classical dynamic scaling is activated [|18| . Extending the argument 
to the quantum dynamics at the T = transition, suggests that the quantum dynamic 



scaling would again be activated |]19[. Further, classical fixed points have also appeared in 
recent studies of the metal-insulator transition PU|, and in the T = onset of spin-glass 
order in metallic systems [|TT|; whether, in these cases, the dynamic scaling is activated or 
conventional, due to the presence of itinerant fermions, is an interesting question, worthy of 
future investigation. 

One of us (T.S) would like to thank S.N. Majumdar for a collaboration on some re- 
lated earlier work |16[. We also thank Y. Gefen, N. Read, and T.F. Rosenbaum for useful 
discussions. This research was supported by NSF Grant No DMR 96-23181. 
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FIG. 1. Phase diagram of the dilute Ising model in a transverse field (h) at zero temperature. 
The dilution probability is p. The multicritical point M is at p = pc, h = Hm- The quantum 
transition along the vertical phase boundary (h < Hm, p = Pc) is controlled by the classical 
percolation fixed point at p = pc, h = 0; quantum effects (due to a non-zero h) are dangerously 
irrelevant, and lead to activated dynamic scaling near the h < Hm, p = Pc line. 
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FIG. 2. Finite temperature phase diagram for h < Km- The dashed Unes (T ~ 1/ ln(l/[p— pd)) 
represent crossovers from the high T regime, characterized by spin fluctuations on the criti- 
cal infinite cluster, to the low T regimes. The solid line for p < pc is the phase transition 
(T = Tc ^ l/ln(l/(pc — p)) where long range order is destroyed. 
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